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We consider a mixture of a superfluid Fermi gas of ultracold atoms and a Bose- Einstein condensate 
■ of molecules possessing a continuous (7(1) (relative phase) symmetry. We study the effects that a 

' spatially random photo-associative-dissociative symmetry breaking coupling of the systems. Such 

coupling allows to control the relative phase between a superfluid order parameter of the Fermi 
j^JTJ' system and the condensate wavefunction of molecules for temperatures below the BCS critical 

r-{ ' temperature. The presented mechanism of phase control belongs to the general class of disorder- 

, induced order phenomena that rely on breaking of continuous symmetry. 
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Of all the studies of the physics of ultracold atoms, 
fermionic systems and in particular Fermi supcrfluids 
have recently received a lot of attention. The most fa- 
mous example in this context is the BEC-BCS crossover. 
In this example, fermions are transformed between the 
superconducting state in which they form Bardeen- 
Cooper-Schriefer (BCS) pairs and the Bose-Einstein Con- 
densate (BEC) made of diatomic molecules of the same 
fermions (for reviews see [H, SQ). Another field of par- 
ticularly high activity has recently been the study of dis- 
order in ultracold gases. Since its proposal [H, several 
experimental and theoretical groups have focused on this 
subject 0, Q. Very recently, these studies culminated 
in the seminal experimental observation of Anderson lo- 
calization of matter waves in a BEC with a disordered 
potential @, S] and in signatures of the Bose glass state 
of an ultracold gas in an optical lattice . 

We have proposed an additional type of disorder ef- 
fects that can be realized with ultracold atoms, called 
disorder induced order [loj ]. Generally speaking, the ef- 
fect can be observed in systems possessing a continu- 
ous symmetry which, in accordance with the Mermin- 
Wagner-Hohenber g th eorem, prevents long range order 
in low dimensions [lH . Comparing models without and 
with continuous symmetries, we therefore observe that 
the latter have lower critical temperatures in any dimen- 
sion. Adding disorder that breaks the continuous sym- 
metry leads to competing mechanisms in such systems: 
naturally, disorder acts against ordering whereas symme- 
try breaking increases the tendency to order. We have 
studied this effect for a classical XY model in a random 
X oriented magnetic field and have shown rigorously the 
appearance of spontaneous magnetization in the Y direc- 
tion at T = 0. Also, we have presented strong evidence 
for the survival of magnetization at T > in the limit of 
small disorder. These ideas were then applied to quan- 
tum systems: a two component Bose gas in an optical 



lattice [10( and a two component BEC with random Ra- 
man coupling [l2| . It is worth noticing that similar effects 
can be found in other areas of physics (see for instance 
[lil] . for classical mean field modelssee [HI)- Also, the ef- 
fect is robust enough to appear not only in random, but 
also in oscillating symmetry breaking perturbations. 

The main goal of the present letter is to introduce 
these new ideas to the field of superfluid Fermi gases 
that posses a phase U(l) symmetry. We show that a 
Fermi superfluid coupled to a molecular BEC via a ran- 
dom, symmetry breaking photo-associating-dissociating 
coupling undergoes relative phase ordering, so that the 
phase of the order parameter can be efficiently controlled 
by the phases of the coupling. Within a BCS-likc theory 
[1 5I | , we show rigorously that for small disorder the effect 
occurs practically for all temperatures below the super- 
fluid transition temperature. 

We consider a mixture of fermions in two different in- 
ternal states interacting via attractive zero-range poten- 
tial in a volume V in 3D with a Hamiltonian 



H F = dr 



$3 



. (1) 



where the chemical potential \x fixes the average density 
n, g = ^-Kh 2 \a\/m (with s-wave scattering length a) de- 
termines strength of the interactions and the ipj stand 
for fermionic field operators. Moreover, we assume the 
presence of a Bose-Einstein condensate (BEC) of molec- 
ular dimers consisting of the two fermionic species and a 
(weak) coupling that transforms them into fcrmion pairs 
and vice versa. Experimentally, we can sweep a mix- 
ture of fermions with two different internal states over 
a Feshbach resonance, which leaves us with a BEC of 
diatomic molecules and unbound fermions. Then, we 
approach a second Feshbach resonance which turns the 
previously unbound fermions into BCS pairs without af- 
fecting the molecular BEC. The coupling between the 
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molecules and the fermions can be realized through pho- 
toassociation and photodissociation, respectively. Tak- 
ing the limit of a large BEC, we do not need to con- 
sider the dynamics of the BEC because the effect of 
the weak coupling with the fermions on BEC is neg- 
ligible. Following these considerations, the Hamilto- 
nian ([I]) has to be supplemented by one term only: 
the coupling between the fermions and the BEC which 

we approximate as / dr 7*(r)^?^i^2 + l{ Y )i } \Hi } d ~ 
y/nd J dr 7*(r)^i^ 2 + 7( r )VWi > where the bosonic 

filed operator ipd for molecules is substituted by a real- 
valued condensate wavefunction which for a homoge- 
neous case considered here is squared root of the density 
of dimers Jri^. The full Hamiltonian therefore reads 



H = H F 



dr 



r*(r)^2+r(r)$#l, (2) 



with r(r) = r(r)e" 4Vr = ^/n^j(r). We assume that the 
transfer process is realized such that T(r) is real, changes 
randomly in space and is constant in time, J c?rr(r) = 
and tpr is a real number. We will show that for pp = 
the relative phase between the condensate wavefunction 
of molecules and the paring function of the superfiuid 
fermions is fixed to tt/2 (or —tt/2). Then we show that 
we can control the relative phase and fix it to any value 
by changing a control parameter pp. 

For r = and in the weak coupling limit, i.e. g — > 0, 
we deal with a Fermi system which for T less than the 
critical temperature T c = 8e ? 'e~ 2 7r~ 1 lVe _7r / 2fcF l a (where 
k B T F = e F = h 2 k 2 F /2m = h 2 (3TT 2 n) 2 / 3 /2m and 77 = 
0.5772) reveals a transition to a superfiuid phase (BCS 
state) which is indicated by a non- vanishing pairing func- 
tion (order parameter) A = g(ip2'>pi} [3 El- I* 1 ^ ne S en ~ 
eral case (i.e. including non-zero T) the pairing function 
is given in terms of solutions of Bogoliubov-de Genncs 
(BdG) equations, i.e. A = g J2 n — 2f(E n )], where 
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= E„ 



with a Hartree-Fock term W 



(3) 



-5E„(K| 2 /(^n) + K| 2 [l- f(E n )}) and f(E n ) = 
( e s n /fe B x + 1 )- 1 

If the transfer process is absent (i.e. T = 0) the system 
(|2"|) is invariant under global gauge transformation, i.e. 
ijjj — > e ltp / 2 ^j, which implies that if {u n , v n } are solutions 
of the BdG equations for A, then {e lv ^ 2 u n , e~ ltp / 2 v n } are 
the solutions corresponding to e %v A. This continuous 
symmetry is broken when the transfer process is turned 
on, as can be seen through ([2]). Then the phase of the 
pairing function becomes relevant because it is the rela- 
tive phase with respect to the (real-valued) condensate 
wavefunction of dimers. 



We begin with an analysis of the pp = case, i.e. 
r = T is real. Let us assume that for F = and for 
some temperature T we have a non-zero pairing function 
which is chosen to be real and positive, A > 0. When 
we turn on V but |r(r)| -c Ao one may expect that it 
results in a new pairing function where A(r) « Aoe zlp ( r \ 
That is, any non-zero T has a dramatic effect on the phase 
because without the transfer process the system is degen- 
erated with respect to the choice of (p. On the other hand 
an infinitesimally small T is not able to change |A(r)| be- 
cause this would cost energy. Moreover, we may expect 
that p(r) oscillates around some average value ipo with 
small amplitude because we assume that T(r) fluctuates 
around zero with infinitesimally small variance. Under 
these assumptions we can observe that \(po\ = tt/2. In- 
deed, let us neglect the Hartree-Fock term W (which is 
not essential for Fermi superfluidity) and calculate the 
difference of the thermodynamic potentials for superfiuid 
and normal phase 
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where H\ = J dr f (V'iV'2 + Vw!) ~ g^li^l^ip 
According to the assumptions, for g' close to g, |A(r)| is 
constant and cos</?(r) « cospo — sin p dp(r). Then 



dr 



,2 , Vft 



A| 2 + — r|A|cos<^ « 
9 

A\ 2 V + sm(p ^^ [ drf6<p, 
9 J 



(5) 



and for J drTdtp < the thermodynamic potential is 
minimized when po = tt/2. With the transformation 
dip — ► —6(f and po —tt/2 we obtain another solution 
what reflects a symmetry of the system. That is, for a 
real T, if A(r) corresponds to solutions of ([3|) then so- 
lutions of complex conjugate BdG equations results in a 
new pairing function equal A*(r). In experiments the 
sign of po w iH depend on the preparation and is deter- 
mined by spontaneous breaking of the p — > — p symme- 
try 

Having determined po we would like to estimate fluc- 
tuations of the phase of the pairing function 6p(r). To 
this end let us employ Ginzburg-Landau (GL) approach 
fljl ]. Adapting the Gorkov's derivation of the GL equa- 
tion [3, [It], EH (with the standard regularization of the 
bare interaction g for the case of cold atomic gases) to 
our problem one obtains 



V 2 A 



-v 2 r 



48tt 2 



7C(3)Z 2 \mkFg 
48tt 2 T c - T A 6m 2 



2tt 2 K 2 



T c -T 



r 



Kim t c 



h 4 k 2 F 



A + r| 2 (A + r), (6) 
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where l c = H 2 kF /mkBT c . Equation (J6]) is valid for T c — 
T<T C and for T(r) that changes on a scale much larger 
than l c (e.g. for kp\a\ = 0.5 and n ~ 10 14 cm -3 we 
get l c ~ 4 fim). For |f (r)| much smaller than |A (T)|, 
i.e. the pairing function in the absence of the transfer 
process, we may introduce further approximations that 
allow us to reduce Eq. ([6]) to 

|A |V 2 Mr) = V 2 f(r) 



T 



r(r), 



where |A| w |A | 
expansion ip(r) ~ 



A8tt 2 

~7((3)P C \mk F g ' T t 

(7) 

and we have chosen ipo = tt/2 in the 
ipo + 5<p(r). The solution of {7J reads 



<^(k) 



f(k) 



48tt 2 



|A | 7C(3)/ C 2 |A | Vmfc F5 T« 



27r 2 n 2 T c -T\ T(k) 



(8) 



in the Fourier space. 

Now we switch to a general case of complex T = 
Te~ lVT . It is easy to check that if lAle 1 ^ corresponds 
to solutions of the BdG equations with ipr = then 
|A|e i(v_¥,r) i s related to the solutions for ip r ^ 0. This 
implies that, if for ipr — we are able to fix the relative 
phase between the condensate wavefunction of molecules 
and the pairing function of the superfluid fermions to 
7r/2 (or —tt/2), then changing ipr allows us to fix it to 
(fio = tt/2 — ipr (or </>o = —tt/2 — ipr) and phase control 
emerges. 

Assuming that the transfer process with small |T(r)| 
results in phase fluctuations of A(r) only, we have shown 
that the fluctuations occur around tt/2 — ipr (or —tt/2 — 
ipr) and they are given by Eq. |8|). Now we would like 
to switch to numerical solutions of the BdG equations 
(where, in contrast to the analytical analysis, we do not 
neglect the Hartree-Fock term W) to demonstrate that 
indeed for |T| -C |Ao| the fluctuations are small and the 
predicted phase control is possible. In 3D calculations we 
regularize the coupling constant g in A = ^(^V'l}) i- c - 
9 -* 9eff, according to 
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Jeff 




2tt 2 H 2 



(9) 



where the logarithm term results from the sum over Bo- 
goliubov modes corresponding to energy above the cut-off 
Ec performed in the spirit of the local density approxi- 
mation, see [la for details. For the simulations we choose 



L, 



40/c; 



which for T 
A (T = 0) 



L± = 20kp L , fi = 0.83e F and fc F |a| = 0.4 
= and the cut-off Ec — IOOef leads to 
= 0.036£p and T c = 0.019T F . Using these 



parameters l c ~ lOOfc^ is larger than the system size 
and we are able to explore a regime beyond GL theory. 
We assume real T(r) given by a pseudo-random function 
that changes along z direction only, 



2tt , 
sin | — (9z 



8.8) 



l( ^(13z + 3.6) 




(10) 



Figure 1: Panel (a) shows T(z) given in Eq. (|10[) for r = 
0.01|Ao(0)|. Panel (b) presents the corresponding phase <p(z) 
of the pairing function for T = (black solid curve) and 
T = 0.9T C (red dashed curve). 



In Fig. [T] we show the phase of the pairing function 
tp(z) in the case when r = 0.01 1 A o (0) | and ipr = for 
two different temperatures, i.e. T — and T — 0.9T C . 
One can see that indeed the phase oscillates around tt/2 
with a small amplitude (standard deviation of the or- 
der of 10~ 2 ). The fluctuations of the absolute value of 
A(z) are negligible (standard deviations divided by aver- 
age values are of the order 10~ 4 ). When T approaches T c 
the average |A| decreases and, at some T, becomes much 
smaller than Tq and we enter another regime where the 
transfer term in the Hamiltonian ((2]) starts dominating. 
For a very large Tq we may expect that a real-valued 
A, which fluctuates out of the phase of T(z), minimizes 
the thermodynamic potential. In Fig. [2] we present aver- 
age values and standard deviations for tp and |A| versus 
temperature where one can observe an increase of the 
fluctuations for T — > T c . 

In summary, we have shown how to control the relative 
phase ip between the wavefunction of a molecular conden- 
sate and the pairing function of a mixture of fermions in 
the BCS state. It turns out that a certain class of weak 
couplings which transfer pairs of fermions into molecules 
and vice versa, fix this relative phase. The couplings need 
to vary randomly, or in an oscillatory manner in space; 
they can be realized by optical means, with a desired 
phase and amplitude, which in turn allows for efficient 
control of ip. In this letter we have considered the Fermi 
system in a weak coupling regime but similar behavior 
is expected in the strong regime. In particular, transla- 
tion of our results to the simplified resonant superfluidity 
theory (cf. (3^) is straightforward. Our results hold also 
for < kpa <C 1, where the pairing function becomes a 
condensate wavefunction of tightly bound pairs and we 
deal with the control of relative phase between two Bosc- 
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Figure 2: Panel (a) shows average value of the phase of the 
pairing function (ip) z versus temperature obtained for V as 
in Fig. [T] In panel (b) we present the corresponding average 
value of the modulus of the pairing function (|A|) Z divided by 
Ao(T), i.e. the pairing for the V — case. Solid black curves 
are related to average values, dashed red curves to average 
values ± standard deviation. The figure shows simulations for 
temperatures up to T = 0.99T C where A (T) = 0.16A (0). 



Einstein condensates, analyzed in our earlier publication 
[13]. The problem considered here belongs to a general 
effect of disorder-induced order phenomena, that rely on 
continuous symmetry breaking. 
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